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Abstract 



In this paper the couphng between dark energy and the other com- 
ponents of the cosmological fluid has been studied. Firstly, it will be 
shown that the application of general cosmological equations, deduced by 
the authors in a previous work, to the known data of the Abell cluster 
A586 using the Layzer-Irvine theory gives similar results compared to the 
work of other authors. The aforesaid method present some problems: the 
application of an approximate theory (Layzer-Irvine theory), it has only 
one experimental datum and finally, it gives results that are a bit difficult 
to admit considering a physical reasoning. In order to avoid the above- 
mentioned problems, a way to study the coupling of dark energy with 
other Universe components has been shown. The results obtained have a 
sensible physical behavior. They also fix the required functionality of the 
product WA^A in order to verify the main known properties of the Uni- 
verse's behavior. Finally, these results permit to make predictions about 
a set of different cosmological properties. 
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1 Introduction 

Some experimental data obtained at the end of the last century, have showed 
that our Universe is speeding up its expansion. Recently, it has been found out 
that more accurate experimental data reach the same conclusion [1-4] 

Several solutions, such as the existence of cosmological constant [5,6], mod- 
ified theories of Gravity [7], dynamical dark energy [8-10] and unified dark 
energy-dark matter models [11], in order to explain the experimental data were 
tried . But, generally speaking, all these models had to face several problems, 
being the main one the coincidence problem and the fine tune of model param- 
eters. 

The newest models have started to study new possibilities like interactions 
between different kinds of components of the Universe, mainly interactions be- 
tween dark energy and dark matter [12-20], new equations of state and mass- 
varying cold dark matter [21]. 

A more general model was developed recently by A.C. Balfagon and R. 
Ramirez-Satorras [22]. In this model some of the assumptions taken into ac- 
count by most authors, were relaxed. For instance, not considering that all 
the Universe's components had no interaction pressure between them. Within 
this recent model, the authors also studied an effective equation of state for the 
whole Universe instead of taking individual equations of state for each of its 
components. 

This paper is focused on the study of the interactions between dark energy 
and the other Universe components, a part from the dark matter. Most of the 
models developed until the present time, have studied the dark energy and dark 
matter coupling; some authors have concluded that there is a energy transfer 
from dark matter to dark energy, while others stand for the opposite. Thus this 
is an issue involving much controversy [23-28]. 

One of the main problems of the interacting models is that, until this mo- 
ment, the interaction has to be fixed beforehand, being the so called scaling 
interaction one of the most popular [29-31]. In this paper a scaling coupling 
and also a new form of the interaction that comes directly from our model have 
been studied. The latter taken into account a sensible physical assumption 
which has never been described before. 
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The rest of the paper is divided in six sections; section two exposes the 
general equations [22], and deduction of new ones. The third section deals 
with a scaling solution, similar to the one used for many authors but using a 
more general reading. The Layzer-Irvine theory is applied in order to find the 
parameters of the coupling using the experimental data from the Abell cluster 
A586. 

In the fourth section, a sensible physical condition is applied to the general 
equations deduced in section two; from this a new coupling function of dark 
energy and dark matter is deduced. In the section fifth diff'erent examples, 
deduced from the model of section fourth, have been studied. Finally, in the 
last section, the main conclusions are exposed. 



2 General equations. 

In [22] general equations for the Universe were deduced. The main idea behind 
the model was to take into account that the components of the Universe could 
interact between them. Thus it was important to consider some kind of inter- 
action pressure between these components. With this idea, the application of 
the deduced equations can be done through the main history of the Universe. 
The main equations are: 

Where p' stands for any kind of energy but the one that come from A, which 
can be considered as a constant (the so-called cosmological constant) or as a 
dynamical function. 

Dark energy density pA will be: 

The total energy density p of the Universe, including dark energy, will be 
defined as: 



The total pressure p of the Universe, without assuming that this pressure is 
only the sum of the partial pressure components, will be: 



(4) 

Equation (5) is the effective equation of state (EOS) for the cosmological 
fiuid (sum of all components). Equation (6) is the function for w proposed in 
[22]. 
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p = wp 



(5) 



1 jl + z) 

To get the best set of parameters in Eq.(6), a full factorial design of experi- 
ments (DOE) [22] was implemented using the data reported in [1]. 

a = 4.1, /3 = 3.5, 7 = 1.0, h = 0.697, H{0) = lOO/i (7) 
With the above equations, the following equations can be easily deduced: 

ia+£) .o. 

^-3(^-(^ 

P = Po[ —5— — - — (9) 
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Where po is the present fluid cosmological energy density. 

The dark energy density variation can be deduced from Eq.(l), expressed 
now using the Hubble function H as: 

H^^^/^A (U) 
The most general variation of A will be, using Eq.(ll): 

A = + + Ao (12) 

Where <J>, \1/ and Aq are considered, in general, not constant. It is easy to 
show that the general expression (12) is equivalent to: 

A = aip' + ao (13) 
Consequently the dark energy density can be expressed as: 

pA = ^A = /3ip' + /3o (14) 

In the last equation /3i and /3o are considered as some kind of effective cou- 
pling between dark energy and the others Universe components. The fact that 
this model includes any kind of coupling apart the one between dark energy and 
dark matter, makes it a newer and more global model. 

With the present data, /3o and /3i, have currently the values and 18/7. 

In the next sections some important models will be studied using these gen- 
eral equations. 
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3 Scaling solution. 



In this section, a scaling solution for f3i in Eq.(14) is studied. This kind of 
coupling has been previously studied by other authors [29-31]; however, consid- 
ering it only as a coupling between dark energy and dark matter. In the model 
described in this article, this scaling solution is taken as an effective coupling 
between dark energy and the other components of the cosmological fluid. The 
Pi function will be: 

In order to compare our results with the ones in [32] we will consider the 

present value of the dark energy density parameter D,\ equal to 0.72, using this 
the value for S has to be 18/7 and (3o has to be 0. Substituting Eq.(14) in Eq.(3) 
the following is obtained: 

PA = (16) 

And substituting Eq.(15) in Eq.(16) the dark energy density evolution is 
deduced: 

18 

= 18 + 7(1 + .).^ 

In order to get the value of the constant r] it is considered that in the present 
period: 

p « PA + Pdm (18) 

Where pdm stands for dark matter density 

With Eq.(17) replaced in Eq.(18) the dark matter density is deduced to have 
the form: 

7(1 + z)^ 

= 18 + 7(l + z).^ ^'^^ 

Replacing Eq.(9) in Eq.(19) we get the dark matter density evolution as a 
function of the shift parameter z: 

_ 7po {z + (3r{a + z){l + zr 
^'"~a/33 18 + 7(1 + ^)" ^ ^ 

Following the procedure developed in [32], and applying the Layzer-Irvine 
equation [33] to the Abell Cluster A586, a value for r] can be obtained. 
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3.1 Layzer-Irvine equation 

The Layzcr-Irvinc equation is the equation for the classical virial equilibrium, 
applied for non-relativistic self-gravitating dust-like particles. It will be assumed 
that the mass of the system remains constant over the observation process. All 
the matter of the system is assumed to be mainly dark matter. The local energy 
density for the dark matter will be the sum of the kinetic energy density and 
the self-gravitating energy density: 

Pdrn = Pk+ pdm (21) 

The pk can be computed considering the Newtonian kinetic energy K per 
unit mass as: 

Where the averaged momentum p and mass M are constants, a is the scale 
factor of the Robertson- Walker's metric. The pk will be: 

Pfc = ^^ (23) 
The pw is evaluated through the potential energy per unit mass: 

W = -27rGaVdm j rC{r)dr (24) 

Where pdm is taken as the background dark matter density, given by Eq.(20), 
and ^ is the auto-correlation function. The pw will be: 

Pw = (25) 
The derivation of Eq.(21) gives the Layzer-Irvine equation: 

- [2pK + (1 + z)Cpw] H (26) 



dPdm 



dt 

Where H is the Hubble function, and C is given by: 

^ z + P z + a 1 + z 18-1-7(1-1-0)'' ^ ' 

If the system is stationary then = 0, then from Eq.(26): 

2pK = -{l + z)Cpw (28) 
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3.2 The Abell cluster A586 



The Abell cluster A586 can be considered as an stationary system, where Eq.(28) 
can be applied, with a redshift of z = 0.1708. The kinetic and gravitational 
energy densities can be estimated as in [32] . The quotient value and band error 
are: 

— « -0.76 ± 0.05 (29) 

Pw 

Using these data in Eq.(27) and Eq.(28) it is possible to get a value for ry, 
this value is 3.95 ± 0.24, while in [32] the value was 3.82. These values indicate 
that the coupling, at least in the recent history of the Universe, between dark 
energy and dark matter is inducing a transfer of energy from dark matter to 
dark energy [32]. 

In Fig.l, the dark energy density evolution for z G [—1,4] using Eq.(17) and 
Eq.(9) is shown. 



Pa/Po 

0.8 r 




-10 12 3 4 

Figure 1: — against z 

It seems clear that the dark energy density behaves in a bit strange way, 
why it is maximum for z w 0?. In the next section the coupling function given 
by Eq.(15) will be changed in order to solve this problem. 

In Fig. 2, the dark matter density evolution for z G [—1,4] from Eq.(20) is 
shown. 
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Figure 2: against z 

Furthemore using Eq.(4) it is possible to see if this dark sector coupling it 
is enough to explain all possible recent interactions between the components 
of the Universe. The pressure in Eq.(4) can be considered as a sum of several 
partial pressure added to a generic interaction pressure: 



P = Wp^p„i+Pr+Ph+ Pi 



(30) 



Where Pm is the matter pressure, Pr the radiation pressure, pA the dark en- 
ergy pressure , Pi the interaction pressure and w is given by Eq.(6). Considering 
the equation of state usually taken for the above components we have: 

Pm = 



0.0001 



{i + zY 



(31) 



PK = -PK 



Where p\ is given by Eq.(17). 

The interaction pressure is deduced replacing the above equations in Eq.(30): 



Pi = Po 



z + py {a + z) 



Po- 



1 + z (3-1 
3{a + z) 13 + z 

0.0001 



18 



18 + 7(1 



(32) 
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In Fig.3 the Eq.32 for z G [-1, 4] is shown. 



0.5 r 




-2.5 L 



Figure 3: — against z 

Watching Fig.3, it seems to be clear that, in the recent history of the Universe 
(for z less or equal to 3) there is still some room for other interactions, or that 
the coupling effective function given by Eq.(15) doesn't work properly in order 
to explain the coupling in the dark sector. An alternative way to explain all the 
interactions for the recent history is showed in the next section. 



4 General coupling function. 

In this section, the possibility of finding a coupling function between dark energy 
and the other components of the Universe using a sensible physical condition, 
will be described. 

The procedure consists in considering Eq.(16) but without specifying any 
kind of function for j3i. Then Eq.(16) is expressed as: 

- YTWi' ^ ^^^^ 

Where, obviously, A is the adimensional density parameter n^. 

The EOS for the cosmological fiuid is given by Eq.(30), where w is given by 
the Eq.(8), p is the Eq.(9) and p,„ and Pr are given by Eq.(31). The EOS for 
the dark energy is now taken as: 

PA = waPa^ waAp (34) 
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Replacing the above equations in Eq.(30) the interaction pressure is deduced 
as follows: 



(z + /3)3 (a + z) / 1 + z 13-1 
P'^ a \6[a + z) p + z 

0-0001 ,, , ^4 

- Po — ^ — (1 + z) (35) 

At this point the physical assumption of considering that, for z < 3, the 
interaction pressure between the components of the Universe is meaningless 
gives rise to: 



^ 1 + z _ _ a/330.0001(l + z)4 

"""^ 3(a + z) + z 3(z + /3)3(a + z) ^ ' 

It is worthwhile noticing that for z = Eq.(36) is equal to —0.63, while the 
most probable value for the adimensional density parameter is 0.72 with w < 0. 
In Fig. 4 the evolution of Eq.(36) for z < 3 is shown. 




Figure 4: wt^A against z 



The high degeneracy between wa and pA is remarkable. In the next subsec- 
tions two important cases will be discussed. 



4.1 = —1 

If wa is assumed to be constant, at least for z < 3 and equal to —1, then from 
Eq.(36) the adimensional density parameter fiA(z) can be deduced: 
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np,{z) = A = 



l + z P^l a/330.0001(l + z)^ 



3(a + z) 13 + z 3(z + /3)3(a + z) 

From Eq.(37) it can be obtained that fiA(O) — 0.63. 

At this moment, the pA is given from Eq.(33) and Eq.(37): 

PA = Apa — (a + z) 



(37) 



(38) 



Fig. 5 shows the evolution of pA- Comparing this evolution with Fig.l it is 
evident that the behavior of Eq.(38) is more physically sensible. 



Pa/Po 
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Figure 5: pa/ po against z 



From the equations (33) and (37) /?i is easily deduced. In Fig. 6 the evolution 
of /3i for z < 3 is revealed. 
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Figure 6: /3i against z 



4.2 Variable scaling function. 

From Eq.(36), and taking r2A(0) = 0.72, the wa deduced is —0.88. Taking wa 
as a constant, at least for ^ < 3, the equations for JIa and for pA are very similar 
to equations (37) and (38). 

In this case it is possible to match Eq.(17) with Eq(33) taking for the A 
function the one deduced from Eq.(36) with wa = —0.88. Both equations give 
rise to the same value for i7A(0). For the match to be possible it is necessary to 
consider the variable r\ not constant, given a variable scaling function: 



Evaluating Eq(39) for the AbeU cluster A586 redshift z = 0.1708, a value 
can be deduce for 77 equal to 1.48; as it can be observed, this is a very different 
value from the one deduce in [32] and from equations (28) and (29). This value 
for 77 indicates that there is a transfer from dark energy to the other components 
of the universe. 

In Fig. 7, the evolution of t;, Eq.(39) for z e [—0.1,3] is shown. It should be 
pointed out that the variable 77 is almost constant and less than 3 (the limit 
value that implies a transfer from dark energy to the other components of the 
Universe < 3 or vice versa 77 > 3, thus the coupling between dark energy 
and the rest of cosmological fluid seem to be inducing a transfer of dark energy 
towards the other components, at least for z < 3. 



r\ = 




(39) 



ln(l + z) 
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5 Conclusions. 



In this paper, the apphcation of a general cosmological model developed by 
the authors previously [22] has been taken into account to study the coupling 
between dark energy and the other components of the Universe. 

Initially, an arbitrary coupling function of the scaling type was assumed. The 
method used, in order to find out the value of the scaling constant rj, was the one 
developed in [32] using the Abell chister A586 data and the Layzer-Irvine theory. 
The results obtained were very similar to the ones deduced in [32]. The value 
for T] was 3.95 which implies that the coupling is transferring energy to the dark 
energy from the other Universe components. Using such coupling function, it is 
clear that there is still some room for other couplings, and that the evolution of 
the dark energy density has a behavior a bit unusual: it grows up to a maximum 
near z — and then goes down. Generally speaking this method has a lack of 
experimental data and suffers from the theory applied: the arbitrary coupling 
function used and the classical Layzer-Irvine theory applied. 

In order to avoid the problems of the above-mentioned method, in the second 
part of the paper the application of a sensible physical condition was consid- 
ered. The assumption adopted was considering that for z < 3 the interaction 
pressure between the cosmological components was negligible and from this it 
was deduced the necessary coupling function in order to fulfill this condition. 

The function deduced from the above condition was the product of wa and 
ft A, which shows the high degeneracy for wa and for Qa- This function has 
some relevant features that made it interesting for the research in dark energy: 

• The value of the function for 2; = is wa^^a = —0.63, while the best 
experimental value for SIa is 0.72. 

• It gives the general form of the product wa^a in order to fulfill two main 
conditions: the accelerated expansion of the Universe and the negligible 
value of the present interaction pressure. 

If WA is assumed to be constant, at least for 2; < 3, then two relevant cases 
are easily studied: 

• Taking wa = —1, it is predicted that riA(O) = 0.63. It is also deduced a 
new coupling function, j3i , between dark energy and the other components 

of the Universe, and the function (f^A = A). From these two functions 
a more sensible function for the dark energy density is deduced. 

• Taking riA(O) = 0.72, the value predicted for wa is —0.88. In this case 
the ^1, Qa and dark energy density functions are similar to the above 
case with wa = — 1. However, now it is possible to deduce a new scaling 
coupling function with rj variable. In this case it is shown that the coupling 
is transferring dark energy to the other components of the Universe, at 
least for the z interval where the assumption that was made is applicable. 
Furthermore, rj varies very little for z € [0,3]. 
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